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Abstract 

The RHIC data on charm production are compared with the fcr-factorization 
approach predictions, both standard NLO QCD and FONLL. The calculated re- 
sults underestimate the STAR Collaboration data. The role of possible nuclear 
effects is discussed. 
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PACS numbers: 25.75.Dw, 13.87.Ce, 24.85.+p 

The investigation of heavy quark production in high energy collisions is an impor- 
tant method for studing the quark-gluon structure of hadrons and the possible nuclear 
effects at early stages of secondary production. The description of hard interactions 
in hadron collisions within the framework of QCD is possible only with the help of 
some phenomenological assumptions which reduce the hadron-hadron interaction to the 
parton-parton one via the formalism of the hadron structure functions. The cross sec- 
tions of hard processes in hadron-hadron interactions can be written as the convolutions 
of squared matrix elements of the subprocess calculated within the framework of QCD 
^ ■ with the parton distributions in the colliding hadrons. 

The most popular phenomenological approach is the NLO QCD collinear approxima- 
tion [1-4], where the cross sections of QCD subprocesses are calculated in the Next-to- 
Leading Order (NLO) of a s series. The Fixed Order plus Next-to-Leading-Log (FONLL) 
[5] also resumes large perturbative terms proportional to a™ log fc (pr/m) with k = n, n — 1, 
where m is the heavy quark mass. In these calculations all particles involved are assumed 
to be on mass shell, carrying only longitudinal momenta, and the cross section is aver- 
aged over two transverse polarizations of the incident gluons. The virtualities q 2 of the 
initial partons are taken into account only through their structure functions. 

The standard QCD expression for heavy quark production cross section in a hadron 1 
- hadron 2 collision has the form 

a i2^QQ = I I dxb . G a/1 {x a , fi 2 F ) ■ G b/2 {x b , /i 2 p) ■ a ab -* Q Q(s,m Q ,fi 2 B ) , (1) 
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where \ip is the QCD factorization scale, x a0 = Atjiq/s, and x b0 = Aijiq/ (sx a ). Here 
G a /i{x a , A*f) and Gb/ 2 (xb, /i F ) are the structure functions of partons a and 6 inside 
hadrons 1 and 2 respectively, and 

a ab ^(s, m Q , fi) = aW R ) ^ + (2) 

is the cross section of the subprocess ab — > QQ as given by standard QCD as a sum of 
LO and NLO contributions [1, 2, 4]. These contributions depend on the parton center-of- 
mass energy s = (p a + Pb) 2 = x a XbS, the mass of the produced heavy quark mq (actually 
they only depend on p = Amg/s), and the QCD renormalization scale /i R . 

The possibility to incorporate the incident parton transverse momenta is refered to 
as the /cy-factorization approach [6-9], or the theory of semihard interactions [10-18]. 
Here the Feynman diagrams are calculated taking into account the virtualities and all 
possible polarizations of the incident partons. In the small x domain there are no ground 
to neglect the transverse momenta of gluons, qy? and q 2 x, when compared to the quark 
mass and transverse momenta, Pit- Moreover, at very high energies and very high p^ 
the main contribution to the cross sections comes from the region of qyr ~ Pit or of 
q 2T ~ p 1T ( S ee [16] for details). The QCD matrix elements of the partonic subprocesses 
are rather complicated in such an approach. We have calculated them in LO. On the 
other hand, the multiple emission of soft gluons is included here. All details of the 
calculations in the /^-factorization approach, presented below, can be found in [18]. In 
our calculations the charm quark mass was taken as m c = 1.4 GeV and we used QCD 
scales /i R = /i F = m T , = m 2 + p\. 

We will firstly consider the cross sections of charm production in pp collisions and 
then we will briefly discuss the situation with nuclear effects. 

The existing data on total cross section of charm production at high energies are 
presented in Fig. 1. One can see the difference in results by PHENIX and by STAR 
Collaborations. This difference is discussed in details in [25, 26]. 

All experimental points except those by STAR Collaboration and the cosmic ray ones 
are in reasonable agreement with NLO QCD calculations (dash-dotted curve), where 
the GRV95 parton distributions [27] were used, and which are compatible with more 
modern analysis (see discussion in [28]). The /^-factorization approach result (solid 
curve) underestimates the data at comparatively low energies because it contains only 
gluon-gluon fusion contribution and is close to the NLO QCD collinear approximation 
curve at higher energies. 

STAR Collaboration obtains also the p^-distributions of D-mesons produced in d+Au 
collisions and scaled to pp collisions at y/s = 200 GeV. These data, taken from [29] are 
presented in Fig. 2. together with theoretical calculations. The NLO QCD result shown 
by dashed curve is in evident disagreement with the data on the same level as in Fig. 1. 
The upper curve of FONLL calculations for charmed quark production taken from [29] 
(dash-dotted curve) also underestimates the data. The result of the /c T -factorization 
approach (solid curve) has reasonable slope in p^-dependence but also underestimates 
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Figure 1: Total cross section of charm production in pp and pp collisions at p > 400 
GeV/c. Experimental fixed target data are taken from [19] and [20], UA2 Collaboration 
data from [21], all they are shown by points. Squares correspond to PHENIX data [22, 
23] and triangle to STAR data [24] . The data of cosmic rays taken from [25] are shown by 
diamonds. The solid curve shows the result of /c^-factorization approach, dashed curve 
corresponds to NLO QCD with only gluon-gluon fusion contribution, and dash-dotted 
curve to total NLO QCD result. 

the data. The level of disagreement in this case is, however, smaller than in the case 
of total cross sections. This should be connected with the fact that a rather large 
contribution to the total cross section comes from low-pr region (p T < 1 GeV/c) where 
both experimental and theoretical uncertainties are rather large. 

It is necessary to note that in Fig. 2 we compare the experimental points for D-meson 
production with theoretical curves for charmed quark distributions. Contrary to the 
mention in [29], together with fragmentation processes where the momentum of D-meson 
is smaller than the momentum of c-quark, there exist recombination processes where the 
momentum of D-meson is larger than the momentum of c-quark. The existence of 
recombination processes in charm production seems to be evident from the experimental 
data on the asymmetry in yields of the so-called favoured and unfavoured D-mesons, 
see discussions in [30-33]. As a matter of fact the produced heavy and light quarks 
have very different transverse momenta but the difference in the components of their 
velocities can be not so large. Possibly, the fragmentation and recombination processes 
in charm quark hadronization balance each other in the processes with not very high pt, 
e.g. the calculated Feynman-x distributions of produced charm quarks in np collisions 
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Figure 2: STAR Collaboration data for p^-distributions of D-mesons produced in d + 
Au collisions and scaled to pp collisions at ^/s = 200 GeV together with calculations 
in fc^- factorization approach (solid curve), NLO QCD (dashed curve), and FONLL [29] 
(dash-dotted curve). 

are in good agreement with the experimental distributions of produced D-mesons (see 
Fig. 5 in [34]). 

The total cross section of charm production was not measured at Tevatron-collider 
energies. However there exist data onpr-distributions of D-mesons at these energies [35]. 
They are presented in Fig. 3 where it is shown that they are in good agreement with NLO 
QCD calculations for charm quarks (dashed curve) based on GRV95 parton distributions. 
The solid curve corresponds to the /c^-factorization approach and it slightly overestimates 
the data, but the agreement should become better in the future when the contribution 
of charmed antibaryons will be added to the experimental data. The results of FONLL 
calculations [29] with fragmentation functions for D-mesons production presented in [35], 
slightly underestimate the yields of D-mesons. Thus all QCD approaches (see also [36]), 
except of the extremal ones, e.g. the one presented in [37], are in reasonable agreement 
with experimental data at y^s = 1.96 TeV. 

Let us shortly discuss the nuclear effects in charm production at RHIC. First of all 
we can say that EMC-effect, i.e. the nuclear deformation of parton distributions [38] 
should increase the total cross sections for charm production, calculated in the NLO 
QCD linear approximation at RHIC energy by 5-10 % [39] in comparison with a linear 
A-dependence. 

The data of PHENIX and STAR Collaborations scaled to binary interactions are 
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Figure 3: The cross sections for D-meson production in pp collisions at y^s = 1.96 TeV 
[35] with \yi\ < 1 together with calculations in /cy-factorization approach (solid curve) 
and NLO QCD (dashed curve). 

presented together in [26]. One can see the absence of nuclear effects from pp to central 
Au-Au collisions from PHENIX data, and from d-Au to central Au-Au collisions from 
STAR data, both on the level of 20% accuracy. Concerning the PHENIX data presented 
in [25] , there is no visible dependence of the charm production cross section on the number 
of collision participant nucleons, N part , in Au-Au collisions at different centralities. 

There exists a factor ~ 2-3 discrepancy [26] in the total cross sections of charm pro- 
duction obtained by PHENIX and STAR Collaborations. Surely, the explanation of this 
discrepancy is completely connected to some experimental problem. However, meanwhile 
experiments don't clarify this point, two possibilities can be imagined: 
i) if we trust PHENIX data, the NLO QCD reasonably describes all experimental data 
except of two cosmic ray points presented in Fig. 1, and then nuclear effects in total 
cross sections of charm production are small; 

ii) if, on the contrary, we trust STAR data, nuclear effects increase the total cross 
sections of charm production about 4-5 times. If the N part dependence of these nuclear 
effects saturates very fast, even the cosmic ray data could be included in the theoretical 
description. The origin of so large nuclear effects can be connected with large non- 
perturbative contributions in high density states (which can be larger than perturbative 
contributions), e.g. string fusion [40], percolation [41], or colour glass condensate effects 
[42, 43] in the interactions with nuclei. 
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In these approaches, above some scale rj c [41] given by the critical percolation string 
density, the strong colour field inside the cluster formed by the overlapping strings pro- 
duces QQ pairs via the Schwinger mechanism as a single string produces light qq pairs. In 
the same way, in the colour glass condensate approach [42, 43] the significant scale is the 
saturation momentum Q s which grows with energy and nuclear size. When Q s > m T (cc) 
, the classical colour field is strong and produce pairs QQ. The production pattern for 
heavy quarks becomes similar to that of the light quarks and an overall enhancement of 
heavy quark production cross section is thus expected. 

Summary : 

1. The /cT-factorization approach predictions are in reasonable agreement with NLO 
QCD for the total cross section of charm production at RHIC energies. We obtain a 
reasonable description of the PHENIX data and we are in contradiction with STAR 
data. 

2. It seems that the main part of our disagreement with STAR data comes from 
low-pr region (p T < 1 GeV/c) where both experimental and theoretical uncertainties 
are rather large. 

3. The predicted pr-distribution of the produced charm at high pt in the kr- 
factorization approach is higher than NLO QCD and FONLL predictions. The kr- 
factorization approach only slightly underestimates the experimental STAR data. 

4. PHENIX data are compatible with the absence of any nuclear effects in charm 
production, whereas STAR data need rather strong nuclear effects at RHIC energies. 
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for useful discussions. This work was supported by Ministerio de Educacion y Cien- 
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